Nonlinear modes and symmetry breaking in rotating double-well potentials 
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We study modes trapped in a rotating ring carrying the self-focusing (SF) or defocusing (SDF) 
cubic nonlinearity and double- well potential cos 2 0, where is the angular coordinate. The model, 
based on the nonlinear Schrodinger (NLS) equation in the rotating reference frame, describes the 
light propagation in a twisted pipe waveguide, as well as in other optical settings, and also a Bose- 
Einstein condensate (BEC) trapped in a torus and dragged by the rotating potential. In the SF 
and SDF regimes, five and four trapped modes of different symmetries are found, respectively. The 
shapes and stability of the modes, and transitions between them are studied in the first rotational 
Brillouin zone. In the SF regime, two symmetry-breaking transitions are found, of subcritical and 
supercritical types. In the SDF regime, an antisymmetry-breaking transition occurs. Ground-states 
are identified in both the SF and SDF systems. 

^ ■ PACS numbers: 42.65.Tg; 03.75.Lm; 47.20.Ky; 05.45.Yv 



INTRODUCTION 



The concept of the spontaneous symmetry breaking (SSB) in nonlinear systems was introduced in Ref. Its 
significance has been later recognized in various physical settings, including numerous ones originating in nonlinear 
optics [2-[E|, Bose-Einstein condensates (BECs) |6j|-[I2> and degenerate fermionic gases [13]. A general analysis of 
the SSB phenomenology was developed too [14[, which is closely related to the theory of bifurcations in nonlinear 
systems jl5j . 

Fundamental manifestations of the SSB occur in nonlinear systems based on symmetric double-well potentials 
(DWPs) or dual-core configurations. A paradigmatic example of the latter in nonlinear optics is the twin-core 
nonlinear fiber, which may serve as a basis for the power-controlled optical switching |2]. DWP settings in optics were 
analyzed theoretically and implemented experimentally in photorefractive crystals [5|. In the realm of matter waves, 
main effects predicted in DWPs are Josephson oscillations Q, the asymmetric self-trapping of localized modes [H, 
qq ■ and similar effects in binary mixtures 0, [To[ • Both the Josephson and self-trapping regimes were implemented in the 
atomic condensate with contact repulsive interactions [ll|. The SSB was also analyzed in one- and two-dimensional 
CS| . (ID and 2D) models of BEC trapped in dual-core configurations [l2| . 

Another dynamical setting which has produced a number of interesting effects in media with the intrinsic nonlinear- 
ity, especially in BEC, is provided by rotating potentials. It is well known that stirring the self-repulsive condensate 
typically leads to the formation of vortex lattices [I?}, although it was found experimentally [18] and demonstrated 
theoretically [19| that giant vortices, rather than lattices, may also be formed under special conditions (when the 
centrifugal force nearly compensates the trapping harmonic-oscillator potential). On the other hand, the rotation 
of self- attractive condensates gives rise to several varieties of stable localized modes, such as vortices, "crescents" 
(mixed- vorticity states), and the so-called center-of-mass modes (quasi-solitons) [20]. Further development of this 
topic was achieved by the consideration of rotating lattice potentials, which can be implemented experimentally as 
an optical lattice induced in BEC by a broad laser beam transmitted through a revolving sieve [2l|, or using twisted 
photonic-crystal fibers in optics [22j. In these systems, quantum BEC states and vortex lattices have been studied 
[23j, as well as solitons and solitary vortices depinning from the lattice when its rotation velocity exceeds a critical 
value 0. 

A specific implementation of the latter settings is provided by the quasi- ID lattice, or a single quasi- ID potential well, 
revolving about its center in the 2D geometry [25|. In particular, the rotation makes it possible to create fundamental 
and vortical soliton in the self-repulsive medium, where, obviously, nonrotating quasi- ID potentials cannot maintain 
bright solitons [26|. Furthermore, the rotation of a DWP gives rise to azimuthal Bloch bands [27l |28|. 
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FIG. 1: (Color online) The pipe waveguide with the intrinsic potential, twisted at rate uo. 



As mentioned above, the static DWP and its limit form reducing to dual-core systems are fundamental settings 
for the onset of the SSB [2j]-[l2j|. A natural problem, which is the subject of the present work, is the SSB and 
related phenomenology, i.e., the existence and stability of symmetric, antisymmetric, and asymmetric modes, in 
rotating DWPs (recently, a revolving DWP configuration was considered in a different context in Ref. [30], as a 
stirrer generating vortex lattices). To analyze basic features of the phenomenology, we here concentrate on the one- 
dimensional DWP placed onto a rotating ring. As shown in Fig. [TJ in optics this setting may be realized as a hollow 
pipe waveguide twisted with pitch 2tt/uj, while the azimuthal modulation of the refractive index, that gives rise to the 
effective potential V(6), is written into the material of the pipe. Alternatively, a helical potential structure can be 
created in a straight sheath waveguide by means of optical-induction techniques, using pump waves with the ordinary 
polarization in a photorefractive material (while the probe wave is to be launched in the extraordinary polarization 
|3l[), or the method of the electromagnetically-induced transparency (EIT) (33j |. including its version recently proposed 
for supporting spatial solitons [34|. In the latter case, one can make the pipe out of Y2SiOs crystal doped by Pr 3+ 
(Pr:YSO) ions [35j . In either case of the use of the photorefractive material or EIT, the helical structure may be 
induced by a superposition of a pair of co-propagating vortical pump waves, with equal amplitudes, a small mismatch 
of the propagation constants k\^ (Afc = k\ — &2 <C fci), and opposite vorticities (±5), which will give rise to an 
effective potential profile, 



where z is the propagation distance, while r and are the polar coordinates in the transverse plane. In terms of the 
BEC, a similar setting may be based on ring-shaped (toroidal) traps, which have been created in experiments [36| 
and investigated in various contexts theoretically [37j. In that case, the rotating periodic potential can be added to 
the toroidal trap [2l| , which is equivalent to the consideration of the rotating ring [28|, [38| . 

In this work, we study basic types of trapped modes and their SSB phenomenology in the ID rotating ring, in 
both cases of the self-focusing and self-defocusing (SF and SDF) cubic nonlinear it ies. In Sec. II we formulate the 
model and present analytical results, which predict a boundary between the symmetric and asymmetric modes, the 
analysis being possible for the small- amplitude potential and the rotation rate close to uj = 1/2. Numerical results 
are reported in a systematic form, and are compared to the analytical predictions, in Sees. Ill and IV for the SF and 
SDF nonlinearities, respectively. The paper is concluded by Sec. IV. 



V(0, z) ~ r s cos (Ak • z + 2S9) , 



(1) 



II. THE MODEL AND ANALYTICAL CONSIDERATIONS 



As said above, we consider the limit of a thin helical shell, which implies a fixed value of the radius in Eq. (pQ), 
r = ro, that we normalize to be ro = 1. Taking the harmonic periodic potential in the form of Eq. (pQ) with 5 = 1, 
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V(9, z) = 2 A cos 2 (0 — ujz\ the corresponding scaled nonlinear Schrodinger equation is 



■ 9 , 



(2) 



where a = +1 and —1 refer to SF and SDF nonlinearities, respectively. Then, we rewrite Eq. (j2j) in the helical 
coordinate system, with 6' = — ujz: 
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■iw^j + 2Acos 2 (0')-a\^ 



(3) 



where the solution domain is defined at —it < 6' < +7r. For the narrow toroidal BEC trap with the rotating potential, 
the respective Gross-Pitaevskii equation, written in the co-rotating reference frame (cf. Refs. [27], [HI), differs by 
replacing the propagation distance, z, with time t [29| . 

Stationary modes with real propagation constant — \i (in terms of the BEC, \i is the chemical potential) are sought 
for as ip (#', z) = exp (—ijiz) (/>(#'), with complex function (p {6') obeying equation 



[Mp 



1 a! 2 . _d_ 



2Acos\0')-cr 



Equation (j3]) conserves the total power (norm) of the field and its Hamiltonian (energy), 



P = 



\W)\de f , 



(4) 
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with the asterisk stands for the complex conjugate. 

The periodic boundary conditions, V{9' -\-2tt) = V{6') and ip(9' -\-2tt) = ip(0'), make Eq. ([3]) invariant with respect 
to the boost transformation, which allows one to change the rotation speed from u to u — N with arbitrary integer N: 



ip (#', z] u - N) = ip (6\ z; uj) exp 



iNO' + i ( -N 2 -Nu)z 



(7) 



hence the speed may be restricted to interval < lj < 1. Furthermore, Eq. ((3]) admits an additional invariance, relating 
solutions for opposite signs of the rotation speed: ip(6' , z\uS) = ip*(0 f , — z\ —uj). If combined with shift uj — >• + 1, the 
latter transformation implies that the solutions for the rotation speeds and 1 — uj (with < uj < 1) are mutually 
tantamount, therefore the rotation speed may be eventually restricted to interval 



< uj < 1/2, 



(8) 



which plays the role similar to that of the fist Brillouin zone in solid-state media [39j. 

An analytical approach for small-amplitude modes can be developed if the amplitude of the potential in Eq. ((3]) is 
small too, \A\ < 1, and uj is close to the right edge of zone (|8j), S = 1/2 — uj < 1/2. In the simplest approximation, 
the corresponding stationary solutions are looked for as a combination of the zeroth and first angular harmonics, 



<p (0 r ) = ao + ol\ exp (if)') , 



(9) 



where ao is fixed to be real, while amplitude a\ may be complex. Indeed, setting uj = 1/2, A = 0, and dropping the 
cubic term, it is obvious that Eq. ((9]) is an exact solution to Eq. (jlj. Then, under the aforementioned conditions, the 
substitution of ansatz (j9]) into Eq. (j4]) yields, in the first approximation (taking care of the balance of the zeroth and 
first harmonics in the equation), the following algebraic equations: 



lido = Aao — crag — 2a|ai| 2 ao, 

\ia\ = Sai — cr\ai\ 2 ai — 2aa^a\ + Aa\. 



(10) 



First, setting a\ = or ao = in ansatz (j9j) corresponds, respectively, to the CW (continuous-wave) mode, and the 
uniform vortical one, with \ip (0 f )\ = const: 



a\ — 0, a = cr (A — fi) . 

n i |2 

a 



0, |aiT =a(A-fi + 5) 



(11) 
(12) 



TABLE I: Different species of stationary modes in the case of the SF nonlinearity (a = +1), labeled by input waveforms which 
generate them. 



Inputs Types of modes 

cos 0' Fundamental-harmonic symmetric (FHS) 

1 + sin 0' Fundamental-harmonic asymmetric (FHA) 

sin 2 6' Second-harmonic symmetric (2HS) 

1 + sin 0' Second-harmonic asymmetric (2HA) 

sm.6' Anti-symmetric (AnS) 

b + sintf', < b < 1 Broken-antisymmetry (BAnS) 



(recall that a 2 = 1). In terms of full equation (j3]) with A ^ 0, these two solutions extend into non-uniform ones 
[with 1^(001 7^ const], which remain (and are categorized as) symmetric with respect to point 6' = 0. The respective 
vortical mode, obtained as the extension of solution ([T2]h is investigated in a numerical form in the next section, under 
the name of FHS (fundamental-harmonic symmetric) mode. 

Solutions to Eqs. ([TO]) with a^a\ ^ give rise to the species of asymmetric modes (categorized as FHA, i.e., 
fundamental-harmonic asymmetric mode, in the next section), with 

at=^(A-n + 2S), \a^ = ^{A-^-8). (13) 

According to Eqs. (j5j) and (|9j), the propagation constant is related to total power (|5]) of soliton (fT3]h P = 27r(ao + |ai| 2 ), 
as follows: 

+ (14) 

In particular, in the case of the SF nonlinearity, with a = +1, Eq. (p3J) demonstrates that the asymmetric mode meets 
the Vakhitov-Kolokolov (VK) criterion, dfi/dP < 0, which is a necessary stability condition for modes supported by 
the SF terms [5o|. On the other hand, in the case of the SDF sign of the nonlinearity, a = — 1, the asymmetric mode 
satisfies the "anti-VK" criterion, d\ijdP > 0, which, as argued in Ref. |4lj |. may also play the role of a necessary 
stability condition, in the respective setting. 

Further, Eq. ([T3|) predicts a transition between the symmetric (FHS) and asymmetric (FHA) vortical modes at 
<2q = 0, i.e., at /i = A+25. Then, Eq. ([T4j) yields the location of this boundary in terms of the total power, P m in = 27ra5. 
In view of the above-mentioned equivalence between the modes pertaining to rotation speeds uo = 1/2 ± 5, the latter 
result predicts the coexistence of the FHS and FHA modes at 

P>Pmin = 27r|(J|. (15) 

Comparison of this prediction with numerical findings is presented below. 



III. NUMERICAL RESULTS FOR THE SELF-FOCUSING NONLINEARITY (a = +1) 

A. Symmetric, asymmetric, and antisymmetric modes 

Solutions to stationary equation (|4]) were constructed by means of numerical code "PCSOM" elaborated in Ref. 
[i^ . In the SF case, the use of different input waveforms makes it possible to identify five distinct species of stationary 
modes, as listed in Table HI 

The symmetry, asymmetry and antisymmetry of the modes designated in Table 1 is realized with respect to point 
6' = 0, while the solution is considered, as defined above, in the region of — tt < 0' < +7r. Further, the fundamental or 
second harmonic ( "FH" or "2H" , respectively) in the nomenclature adopted in the table refers to a dominant term in 
the Fourier decomposition of the stationary solution (which is made obvious by their shapes, see below). In particular, 
the stationary patterns of the FHA and 2HA types are generated by the same input in Table 1, 1 + sin^, but in 
non-overlapping regions of the parameter space, (cj, A, P), as shown below. Note also that all the inputs displayed in 
the table are real functions, but the numerical modes found for uj ^ have a complex structure. 

The stability of the stationary solutions has been identified through the numerical computation of eigenvalues for 
infinitesimal perturbation modes. To this end, the perturbed solution was taken in the usual form, 



ijj = e-^ z [(j)(e') + u(0')e iXz + v*(0*)e- iX * z ], 



(16) 
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FIG. 2: (Color online) Examples of stable FHS and FHA (fundamental-harmonic symmetric and asymmetric) modes, in the 
case of the SF nonlinearity, corresponding to the same parameter set, (cj, A, P) — (0.35, 0.5, 1). (a,b) Real and imaginary parts 
of the FHS and FHA stationary solutions, respectively. (c,d) Intensity profiles of the same solutions. The dashed curves in (c) 
and (d), and in similar panels displayed below, depict the potential, V {0') — 2Acos 2 (0'). 
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FIG. 3: (Color online) Examples of stable modes of the 2HS and 2HA (second-harmonic symmetric and asymmetric) types, in 
the case of the SF nonlinearity, found at sets of the parameter values (lj,A,P) = (0.25, 1,0.2) and (0.25, 1,0.3), respectively. 
(a,b) Real and imaginary parts of the 2HS and 2HA stationary solutions, respectively. (c,d) Intensity profiles of the same 
solutions. Unlike the fundamental-harmonic modes shown in Fig. [2j these ones, based on the second harmonic, feature no 
bistability (coexistence of the stable modes). 



where u{0') and v{0') are perturbation eigenmodes, and A is the corresponding eigenfrequency. The substitution of 
expression ([T6]) into Eq. (|3j) and linearization leads to the eigenvalue problem in the following form: 

{^ttr -» + t^) {:)--<:)■ <-> 

where h = — (l/2)d 2 / d6' 2 + 2Acos 2 (0 / ) — 2cr\(j)\ 2 is the respective single-particle Hamiltonian. The underlying solution 
cj) is stable if all the eigenvalues are real. 

In the case of the SF nonlinearity, the numerical solutions reveal the existence of all the species of the modes 
indicated in Table 1, except for BAnS, which exist under the SDF nonlinearity, see below. Typical examples of the 
five species of the stationary modes which are supported by the SDF cubic term are displayed in Fig. [2JHJ The 
analysis demonstrates that all these examples are stable. Moreover, the symmetric and asymmetric modes dominated 
by the fundamental harmonic, FHS and FHA, demonstrate mutual bistability in Fig. [2j as these stable modes coexist 
at common values of the parameters. 

A distinctive feature of the FHS mode, which is evident in Fig. [2j is that maxima and minima of its intensity 
coincide with local maxima and the minima of the harmonic potential, while its FHA counterpart features an intensity 
maximum in one potential well, and a minimum in the other. This structure of the FHS suggests that it may correspond 
to a maximum, rather than minimum, of the energy (see below), but, nevertheless, this mode has a region of stability 
against small perturbations. 

On the other hand, Figs. [3] and |4] demonstrates that both the 2HS and AnS modes have two symmetric intensity 
peaks trapped in the two potential wells (while local minima of the intensity coincide, quite naturally, with potential 
humps), hence these mode have a chance to realize a minimum of the energy (in the best case, it may be the system's 
ground state). The 2HA mode features a similar property in Fig. [3j but with unequal density peaks trapped in the 
two potential minima. 

In the limit of the uniform ring (without the potential, A = 0) the FHS and AnS modes go over into the above- 
mentioned uniform vortical state (fT2]h the 2HS pattern degenerates into the CW state (fTTjh while the FHA mode 
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FIG. 4: (Color online) An example of a stable mode of the AnS (antisymmetric) type, at (cj, A, P) = (0.25, 1, 0.3): (a) real and 
imaginary parts of the stationary solution; (b) its intensity profile. 




FIG. 5: (Color online) (a) The FHS mode, in the case of the SF nonlinearity with (cj,P) = (0.25,0.2), at different values of 
the potential's strength, A. (b) The same for the FHA mode with (oj,P) = (0.5,0.5). (c) The same for the 2HS mode with 
(cj, P) = (0.25, 0.2). (d) The same for the 2HA mode with (cj, P) = (0.4, 0.2). All the modes shown in this figure are stable. 



takes the form of a 27r-periodic cnoidal-wave solution of the nonlinear Schrodinger equation equation (the 2HA pattern 
does not exist at all in the limit of A = 0). The evolution of the shape of the modes with the increase of A is shown 
in F ig. [5] [The evolution figure of AnS mode, which is not shown here, is similar to the 2HS type displayed in Fig. 



B. Existence and stability diagrams for the different modes 



Results of the systematic numerical analysis are summarized in Fig. [6l in the form of diagrams for the existence of 
the stable FHS and FHA modes in the plane of (P, uo) at several fixed values of amplitude A of the harmonic potential. 
In the blank areas, only unstable modes of the FHS type are found [in direct simulations they feature an oscillatory 
instability, which is accounted for by a quartet of complex eigenvalues generated by Eq. ([TT]) ]. while the FHA modes 
are completely stable in their existence region. 

The dashed lines in panels (a) and (b) of Fig. [6] demonstrate that Eq. (|T5j) correctly predicts the bistability 
boundary at small values oi 5 = 1/2 — uo and P, provided that A is small enough too [recall it is the condition under 
which Eq. ([T5]) was derived]. Further, the fact that vertical cuts of the panels in Fig. [6j that can be drawn through 
uo = const, go, with the increase of the total power (P), from the region of the monostability of the FHS mode through 
the FHS/FHA bistability area into the monostability region of the FHA mode (at least, for uo sufficiently close to 
1/2), clearly suggests that the symmetry-breaking bifurcation of the subcritical type, which features the bistability 
[l5|, Q, occurs along this route. This conclusion is confirmed by a direct consideration of the numerical data (not 
shown here in detail). 

As concerns the AnS solitons, in the case of the SF nonlinearity they do not undergo any bifurcation, and are stable 
in the entire region of their existence, which is shown in the (cj,P) plane for fixed values of the lattice's strength, A, 
in Fig. [3 [as mentioned above, at A = the AnS mode coincides with the FHS one, hence the respective existence 
area is the same as the red (bottom) one in Fig. [6ja)]. From here we conclude that the stability region of the AnS 
mode quickly expands with the increase of A, and the FHS/FHA bistability areas in panels (b) and (c) of Fig. [6] 
actually support the tristability, as the AnS mode is stable there too. 

The multist ability between the three species of the modes based on the fundamental harmonic — symmetric, asym- 
metric, and antisymmetric ones — makes it necessary to compare their energies, which are defined as per Eq. (J6j). The 
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FIG. 6: (Color online) Stability diagrams for the symmetric and asymmetric modes dominated by the fundamental harmonic 
in the case of the SF nonlinearity, in the plane of the rotation speed (cj) and total power (P), at several fixed values of the 
potential's strength: (a) A = 0, (b) A — 0.2, (c) A — 0.5, and (d) A = 1. The green (top) and red (bottom) areas designate 
regions where solely the FHS or FHA mode is stable, respectively. The bistability (coexistence of these stable modes) occurs in 
the yellow (middle) region. In the blank area, no stable modes of these types are found (in fact, an unstable FHS mode exists 
in this area). In panels (a) and (b), the dashed lines represent the analytical prediction ([T5]) . P — 2n\l/2 — cj|, for the existence 
boundary of the FHA mode (the prediction is relevant at small values of 1 1/2 — uj\, P, and A). The horizontal dotted line in 
(c) is the cross section of the bistability region, along which energies of the FHS and FHA modes is compared in Fig. IIII Bf a). 
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FIG. 7: (Color online) In the case of the SF nonlinearity, the AnS (antisymmetric) mode exists and is stable in the red area, 
as shown here for (a) A — 0.2 and A — 0.4 (b). 



comparison along the horizontal dotted line, which is drawn in Fig. |6(c)| is presented in Fig. 
shows the following relation between the energies: 



a) The figure clearly 



#FHA < #AnS < #FHS- 



(18) 



For the 2HS and 2HA patterns, which are based on the second angular harmonic, Fig. [9] displays the stability areas 
in the (P, uj) plane, at different values of the lattice's strength, A [cf. Fig. [6], together with the stability area for the 
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1.25, 0.5), in the interval of 0.2 < u < 0.5. (b) The energies of the 2HS, 2HA, 



FIG. 8: (Color online) (a) Ene rgies of the FHS, FHA, and AnS modes, computed according to Eq, 
dotted line drawn in Fig 



u with (P, A) 

and FHA modes, along with the coexisting AnS one, along the horizontal dotted line drawn in Fig. 9(c)| for (P, A) — (0.5, 0.5). 
Panels (a) and (b) pertain to the system with the SF nonlinearity, a = +1. (c) The energies of the AnS, BAnS (broken- 
antisymmetry), and 2HS modes, along the boundary separating the existence regions of the AnS and BAnS solutions in the 
case of the SDF nonlinearity (a = — 1), see Fig. 12(b) 
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FIG. 9: (Color online) The same as in Fig. [6j but for the following types of the solutions. The green (top) area: the stable FHA 
mode; the red (bottom) area: the stable 2HS mode; the blue (intermediate) area: the stable 2HA mode; the yellow (rightmost) 
area: the bistability of the FHA and 2HS modes. 



FHA mode, which may stably coexist with the 2HS pattern. Figure [9] also demonstrates that, as mentioned above, 
the 2HA mode, which exists and is stable in the blue (intermediate) area in Fig. [9j does not exist at A = 0, emerging 
and expanding with the increase of A. Note that, on the contrary to the situation for the patterns dominated by the 
fundamental angular harmonic (FHS and FHA) described above, there is no overlap between the stability areas of 
the symmetric and asymmetric second-harmonic-based modes, 2HS and 2HA, hence (as confirmed by an additional 
consideration of the numerical results) the symmetry-breaking 2HS — )• 2HA transition, following the increase of P at 
amounts to a supercritical bifurcation, which features no bistability [TBI, Q. 
Another peculiarity of the symmetric mode based on the second harmonic (2HS) is that it is stable in two disjoint 
stability areas — the red and yellow ones in Fig. [9] (the bottom and rightmost regions), in the latter one the 2HS 
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(a) (b) 

FIG. 10: (Color online) A typical example of a stable BAnS mode (the one with the broken antisymmetry) in the system with 
the SDF (self-defocusing) nonlinearity, a = — 1, and parameters (lj,A,P) = (0.2, 1,0.5): (a) the real and imaginary parts of 
the stationary solution; (b) its intensity profile, which features unequal peaks in the two potential wells. 



featuring the bistability with the FHA mode, although they are not linked by any bifurcation. On the other hand, 
there is no overlap between the stability areas of the asymmetric modes based on the different angular harmonics, 
FHA and 2HA, therefore these modes may be generated by the same input waveform (1 + sin^, in Table 1), in 
different parameter regions. 

The effective bistability occurring in Fig. [9] makes it necessary to compare the energy of the coexisting stable 
patterns, which is done in Fig. |8(b)| along the horizontal dotted line drawn in Fig. [9fc). The energy of the coexisting 
stable AnS mode is included too. In particular, the short segment corresponding to the 2HS mode appears, above the 
one pertaining to the FHA mode, when the horizontal dotted line in Fig. [9fc) enters the yellow (rightmost) region of 
the 2HS-FHA bistability. Thus, adding the results from Fig. |8(b)| we extend energy relations ([18]) as follows: 

i^FHA Or #2HA < #2HS < #AnS < #FHS- (19) 

The conclusion is that the asymmetric modes, either FHA or 2HA, represent the ground state of the rotating ring 
carrying the harmonic potential and SF nonlinearity. 



IV. NUMERICAL RESULTS FOR THE SELF-DEFOCUSING NONLINEARITY (a = -1) 

In the case of the SDF nonlinearity, symmetric modes of the FHS and 2HS types have been found too, although, on 
the contrary to the system with the SF sign, they do not undergo any bifurcations, hence the FHA and 2HA species do 
not exist in this case. On the other hand, the AnS mode features an antisymmetry-breaking bifurcation, which gives 
rise to a broken- antisymmetry (BAnS) mode (which does not exist under the SF nonlinearity). A typical example of 
the stable BAnS mode is displayed in Fig. [10J In fact, antisymmetry-breaking bifurcations are characteristic to the 
double- well systems with the SDF nonlinearity The stationary solutions of the BAnS type can be numerically 
generated by the initial guess of the form of (b + sinfl'), with < b < 1, see Table 1. 

Proceeding to the systematic description of the modes supported by the SDF nonlinearity acting in the combination 
with the rotating harmonic potential, we note, first, that the numerical analysis demonstrates the existence and 
stability of the mode of the 2HS type at all values of (cj,P, A), therefore this solution is not included in the stability 
diagrams, which are displayed for the FHS mode in Fig. [12j and for the AnS and newly found BAnS ones in Fig. [TTJ 

In the blank area of Fig. [12j the FHS solutions exist too, featuring an instability accounted for by a quartet 
of complex eigenvalues [see Eq. (fT7|) ]. In direct simulations, this instability converts the stationary mode into an 
oscillatory one (not shown here). Further, Fig. [T2l shows that the antisymmetry-breaking transition, which follows 
the increase of the total power, P, is of a supercritical type. Actually, the supercritical character of the antisymmetry- 
breaking bifurcation is a characteristic feature of systems with the SDF sign of the nonlinearity [8]. 

Finally, the energy comparison for the model with the SDF nonlinearity is presented in Fig. 8(c)| along the boundary 



separating the existence (and stability) regions of the AnS and BAnS modes in Fig. \Mbj~ The fi sure shows the 
energies of these modes, which coincide along the boundary, along with the results for the FHS and 2HS symmetric 
modes. From here, it is concluded that the energies are related as follows: 

#2HS < #BAnS = #AnS < #FHS, (20) 



i.e., the 2HS mode represents the ground state in the case of the SDF nonlinearity, cf. Eqs. ([T8]) and (p~9j). 
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0.0 0.1 




FIG. 11: (Color online) The stability diagram for the FHS mode, in the case of the SDF nonlinearity (a = —1), in the (lj, P) 
plane, for A = 1.2 (a) and A = 1.75 (b). The FHS modes are stable and unstable in the red and blank areas, respectively. 




(a) (b) 



FIG. 12: (Color online) The stability diagrams for the AnS (antisymmetric) and BAnS (broken- antisymmetry) modes in the 
case of the SDF nonlinearity for A = 0.5 (a) and A = 1 (b). The red (bottom) and green (top) areas refer to the monostability 
of the AnS and BAnS modes, respectively. The bsitability between these species occurs in the yellow (rightmost) region. 



V. CONCLUSION 

The objective of this work is to study the trapped modes of the symmetric, asymmetric, and antisymmetric types, 
and the symmetry- or antisymmetry-breaking transitions between them, in the rotating ring carrying the DWP 
(double-well potential) and the cubic nonlinearity with the SF or SDF (self-focusing or defocusing) sign. The analysis 
has been performed in the first rotational Brillouin zone. In the SF case, five types of different modes, and their 
stability, have been identified, three dominated by the fundamental angular harmonic (FHS, FHA, and AnS), and two 
others based on the second harmonic, 2HS and 2HA modes. The SSB (spontaneous symmetry-breaking) transition 
between the FHS and FHA modes is of the subcritical type, featuring a conspicuous bistability area, while the 
transition between the between the 2HS and 2HA states is supercritical. There is no overlap between the existence 
regions of the asymmetric modes of the FHA and 2HA types, one of which realizes the ground state of the SF system. 

The SDF model supports four distinct species of the trapped modes, viz., FHS, 2HS, AnS, and, in addition, 
the BAnS (broken-antisymmetry) mode. The latter one appears from the AnS state as a result of supercritical 
antisymmetry-breaking transition. The 2HS mode represents the ground state of the SDF system. 

The present work can be naturally extended in other directions. First, it is possible to consider the rotating potential 
in the form of A cos (n6 r ) with n > 3, while the present analysis corresponds to the DWP with n = 2 [note that Eq. 
(PQ) suggests a possibility to create the helical photonic lattice with n = 2S and S > 1]. This generalization will give 
rise to many new trapped modes; in particular, it may be possible to build one with a small width, AO' <C 27r, hence 
it may be considered as a soliton trapped in the lattice and rotating along with it, cf. Ref. [[24']. It is relevant too 
to consider effects of the two-dimensionality (the entire rotating plane, rather than a narrow ring). Our preliminary 
analysis, based on simulations of the 2D counterpart of Eq. (|3]), demonstrates essentially the same types of trapped 
modes in rotating annuli of a finite radial size, as reported above in the one-dimensional limit. Another interesting 
extension is to replace the linear rotating potential by its nonlinear counterpart, generated by modulation of the local 
nonlinearity along the angular coordinate (as was done in Ref. [43| in the absence of the rotation). Results obtained 
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for one-dimensional modes trapped in the rotating nonlinear potential will be reported elsewhere. 
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